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Abstract 

This second part is devoted to the investigation of global properties of Prescribed 
Mean Curvature (PMC) foliations in cosmological spacetimes with local U{1) x U{1) 
symmetry and matter described by the Vlasov equation. It turns out, that these 
spacetimes admit a global foliation by PMC surfaces, as well, but the techniques to 
achieve this goal are more complex than in the cases considered in part I. 
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1 Introduction 



For the main motivation see the introduction in part I, [[H3[| . Only some remarks, special 
to the present situation remain. 

The structure of this second part is as close as possible to the structure of part one. 
Although the preliminary section 2 of part I has been omitted, the formulas here refer to 
it as well. 

The spacetimes considered here are cosmological spacetimes with two commuting local 
Killing vector fields. This symmetry will be referred to as (local) U{1) x U{1) symmetry. 
In comparison with part I, there are now only two (local) Killing fields instead of three. 



generalizing the plane symmetric case of |H3]. The absence of a third Killing field requires 



a more detailed description of the geometry and a deeper analysis to control the momenta 
of the Vlasov particles. Thus the estimates in this work rely on the simple structure of 
the Vlasov equation. Despite this there seems to exist no crucial obstructions for other 
types of 'well behaved' matter. 
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2 Spacetimes with local U{1) x ^7(1) symmetry 



2.1 The geometry of spacetimes with local U{1) x U{1) symmetry 

Following the analysis in we consider now the globally hyperbolic spacetime {M,g) 
with topology M x S, where S denotes a bundle with base and compact orientable 
fibre F. As usual it is assumed, that the submanifolds {t} x S are Cauchy hypersurfaces 
in M. The coordinates of S are denoted by {x,y'^,y^) = {x,y^), A = 2,3, where (y^) 
denote coordinates on F. As usual, Greek indices range in the interval 0, . . . , 3, lower case 
Latin indices from the middle of the alphabet take values in 1, ... ,3, while those from 
the beginning of the alphabet take the values 0, 1 and upper case ones are confined to the 
values 2,3. 

The covering map R — > defines a pullback of the bundle S with base M, hence we get 
a trivial bundle. If F denotes the universal covering space of F, we get a natural covering 
S = M X F of E with canonical projection p. Now we can associate a spacetime {M,g), 
where M = M x S and g is the pullback of g under the projection id xp. 

The fibres F in the trivial bundle M are assumed to be the orbits of a two dimensional 
translation group G of isometries of g. Hence F is the Euclidean space form E2 and 
F = E2/T, for a discrete subgroup F of G. The compactness and orientability of F 
implies then i*" = T-^, so F can be represented by a two-parameter lattice, S turns out to 
be a torus bundle over the circle and the induced action of G on {M,g) is given by the 
quotient action G/T = U{l)xU{l) with orbits F. As in [ [H3| ] we will call the orbits surfaces 
of symmetry, hypersurfaces in M diffeomorphic to {t} x S, which consist of a union of 
surfaces of symmetry will be called symmetric surfaces and we call {M,g) a spacetime 
with local U{1) x U{1) symmetry. 

The induced action of U{1) x ^7(1) on {M,g) is local, if the bundle is non-trivial. In 
that case we have to deal with non-trivial transformations for the transition x i — > x + 2it 
in by lattice preserving translations and automorphisms GL(2, Z) of G. To represent 
the metric g of the orbits in M at first locally for x £ [0,27r[, define the area radius by 
r := (47r)~-^ Vol(-F), and write the metric as 

(1) 9 = r^~9 , 

with a metric g of unit determinant. Due to symmetry both quantities, r and g do not 
depend on the points of F, and one easily verifies, that the curvature of g vanishes, as 
required. For g there are two remaining degrees of freedom, V and W, and we can use 
them to parametrize the metric: 

, s ~ — ( cosh V sinh V 

^ ^ ^~ \ sinhl/ e-^coshy 

If the bundle is trivial, then this representation of the metric is global, but if the bundle 
is not trivial, then the translation of 27r in induces a transformation of g by an element 
Z of GL(2,Z), 

(3) g{x + 2TT) = Z^g{x)Z , 

which fixes the geometry of the spacetime. 

Given a globally acting symmetry we can specialize to some well known geometries: If 



V = W = then we get the plane symmetric case e = of |lH3| . Setting only V = 
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we get a symmetry called polarized, corresponding to the reflection symmetries i — > 
—y^ or i — > —y^ respectively. If the composition of these reflections is an isometry 
(regardless whether the individual reflections act isometrically), we call this symmetry 
to be of Gowdy-type, since the Gowdy spacetimes are defined by this symmetry and 
the additional requirement that the spacetime is vacuum. Thus, in general V ^ m a 
spacetime with Gowdy-type symmetry, but if V vanishes also we call the Gowdy-type 
symmetry polarized. 

We construct the coordinate system {x'*} mentioned in the beginning of this section locally 
by the following procedure. Consider first an arbitrary symmetric Cauchy surface S in M. 
Then we extend the coordinates (y^) of F to a Gaussian neighbourhood of F in S. Later 
on we will do some rescaling along S^, such that we choose the coordinates, such that the 
metric h of S takes the general form h = A"^ dx + g. Now we embed this structure into 
the spacetime. In general the time coordinate t defines lapse N and shift v = di, thus 
the spacetime metric g is not in block diagonal form. The components goA which prevent 
g being block diagonal can be represented by i^a '■= 9ab ■ 



2.1.1 The 2-1-2 geometry 

Now we want to investigate the geometry of {F, g) in (M, g) more closely. As before let 
S denote a symmetric Cauchy surface in M, foliated by its symmetric surfaces. The unit 
normal of F in 5 C M will be denoted by m and n is the unit normal of S in M, as usual. 
Then there are some canonical geometrical objects induced: 



(a) the two second fundamental forms 

(4a) A(v, w) = gCVyW, m) 
(4b) k{v,w) = g{^VvW,n) , 



for i;, u; G TF with arbitrary extensions to vector fields, in order to make the expres- 
sions well defined. In the sequel it turns out to be convenient, to deal with the trace 
free parts A and k instead, defined byA = A + itrA5 and k = k + ^ir ng respectively. 
In the given coordinates we have explicit formulas: 



><AB = -\rn{gAB) = -\A ^{gAsY 

trA = -la^^migAB) = 'h'^^^^f = -'Mr) = -^A'^^ 

KAB = -^n{gAB) = -\N-^{gAB)' - N'^Av^Xab 
trK = -^n(r) , 

reflecting the definition of the area radius r as a volume measure and the second 
fundamental forms as its rate of change. In the above formulas g is considered as 
intrinsic to F (and we maintain this from now on) and g^^ = r~^g^^ denotes the 
inverse of gAB- Differentiating g directly one gets from its relation to g and the 
formulas above the following representations of the trace free parts of the second 
fundamental forms: 

~Xab = -ym{gAB) A^^ = \r-^ m{~9'''') 

kAB = n{9AB) = n{g^^) 
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(b) the connection in the normal bundle T-^F. This can be represented by a single one 
form r], defined as 

rj{v) = g{^V^n,m) = -g{^Vvm,n) 

(5) 1 v& TF 

= --;^9{[n,m\,v) 

From this formula one can see, that = is equivalent to [n, m] G T^F, thus 
[T-^F, T-^F] C T^F and the theorem of Frobenius tells us, that T^F is an integrable 
distribution of two planes in TM. If T^F is integrable, then we can decompose the 
spacetime into a direct sum {M,g) = {F^,g) + {F,g), with g = g ®g. Sufficient for 
the existence of an integral manifold F-^ of the distribution T-^F is the Gowdy-type 
symmetry i — > —y^, A = 2, 3, since then it follows ry = immediately, because 
V I — > r]{v) is an antisymmetric map. 

Whether or not T-^F is integrable, we can orthogonally split the tensor bundle over M, 
following Kundu ( @ ) • To perform this task we start with the two dimensional Riemannian 
manifold {F,g). With our choice of spacetime coordinates we have Killing fields Ya = i9a- 
Their spacetime components define projection operators 

PA^ ■■= Ya^ , PA = dA 

•= 9^^9^.uPB'' , P^ = i^'^dt + dy^ , 

where v"^ := g^^vs = ^"^i and the metric components giB are zero by our definition of 
the coordinates. 

The projection operators into T^F can now defined as follows: Let the unit normal 
vectors n and m serve as projection operators qa^ and define their duals by the relations 
Qa^q^ ^ = 5^ and PA^q^ ^ = 0. The result is 

qQ = n = N~\dt-v'di), q^=m = A-^di 

q^ = Ndt, q^ = Ah'^dt + Adx 

Now we are ready to define the transversal metric g in T-^F by 

(6) gab = gfiuQa^qb" , 

thus in the given frame g and its inverse are represented by the two dimensional Minkowski 
metric and from now on we consider g as an intrinsic object in the tensor bundle over T-^F. 
In summary we have constructed a complete set of projection operators, projecting or- 
thogonally tensors over {TM,g) into {TF,g) and {T-^F,g), characterized by the relations 

5^ (completeness) , 
PA^q'^ = 

^a'^gV = (orthogonality) , 

with PA = Oa and go -L Qi (these two relations fix the component representations shown 
above) . 

Now we define the convention already used in some expressions above, that wherever 
confusion might arise, we attach a bar to quantities, which will be considered as intrinsic 
to the associated bundle, and a tensor index furnished with such a bar denotes projected 



Pa'p"", = 6^ 

qa^p'^a = 
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components, which also can be considered as intrinsic. For example T^b denotes a tensor 
in the bundle {T-^F, g), but Tab = Qa'^Qb'^Tnu, too, although T might be a spacetime tensor. 
The philosophy lying beyond this notation is, that quantities with a bar attached to them 
or to their indices can be manipulated by the associated metric, while indices without a 
bar always denote component indices corresponding to the bundle the tensor is intrinsic to. 
If we apply this notational convention to the metric itself we get the definitions for g and g 
back, for example we get the identities gab = qa'^Qb'' Qij-u = 9ab and gjB = PA^PB^g^u = Qab 
In particular we have the relations 

defining the Levi-Civita connection in the projected bundles, and the algebraic identity 
Ti: = n + Ti . 

For later use we need the projected components of the Ricci tensor. First note, that the 
normal connection -q on {F, g) reads in our new notation 

= -2 bo,gi]^ = -\p^^[qo,qiT = 2^(^0,51) , 

or cquivalcntly ea6^^ = ^/'^dp'^{qa,qb), making use of the relation = {p"^ f^qa^) , which 
gives Vga^ = —qa'^Ps^^P^ v- '^^'^ ^ab here has its standard meaning as the totally anti- 
symmetric symbol (independent of the frame used, so it is not necessary here to attach 
a bar to it). On T-^F we find dp a = d{gABP^) = 9ab dp^ , so that index manipulations 
with g can be applied, as desired. 

After some calculations we arrive at formulas for the projected components of the Ricci 
tensor: 

(7a) 'Rob = gab - 2r-^VaYbr + 2r-\Yar){Ybr) + l{YagAB){Ybg^^) + ^gablvf 
(7b) 'Rjb = r-\%Ycir^VA) 

(7c) 'Rab = -kgADY^ig^'^YcgBc) - r-\Y"r){Yc9AB) - IrjAfiB 

Contraction of the last equation yields a formula for the Laplacian of r (Note, that the 
Laplacian is really a wave operator in the present situation) : 

(8) Ar = + r-HVcr)(V"r) + r|??p 

Finally we consider other frames of reference and their relation to the frames already in use 
in the projected spaces. At first we concentrate on {TF,g), which is canonically endowed 
with the coordinate vector fields Oa- Alternatively we can introduce the orthonormal 
frame {ca}, defined by ca ■= r~^eA, with 

62= e-^/2coshy/2 ^2 - e^/2sinhy/2 ^3 
63 = -e-^/2gij^l^y/2 52 + e^/2poshy/2 ^3 . 

The dual frame {a"^} is defined by := ra^, with {d-^)ij,{eB)'^ = S^, thus 

^2 =e^/2 cosh V/2 dy'^ + e'^''^ sinh V/2 dy^ 
^3 =e^/2sinhF/2 V + e-^''^coskV/2 dy^ . 
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Finally we introduce the two dimensional quotient manifold B := M/(U{1) x U{1)). Then 
we have isomorphisms TB ~ TM/TF ~ T^F and we can consider the metric g as acting 
on TB, so having constructed the two-dimensional (quotient) spacetime {B,g). If we take 
{t,x) as coordinates of B, then the metric g has coordinate components gab = —{q^)'^ + 
(q^)"^ = {—N"^ + (Au^)'^) dt^ + 2A^v^ dtdx + A^ dx"^, giving an alternative way to describe 
tensor components in (T^F^g), even if this bundle is not integrable. The dual base {cq} 
to {a"' = q"'} in {B,g) has the coordinate components cq = N~^{dt — i^^ dx), ei = A~^ dx- 
The component notation just introduced conflicts with the conventions described above, 
and in the following sections we are forced to make clear, which conventions we will follow. 

2.1.2 The 3+1 geometry 

To describe the 3-geometry of an arbitrary symmetric Cauchy surface S, we construct 
the coordinate system {x', y^) explicitly. On the surface of symmetry F C S we have 
already coordinates, such that the scaled metric g takes the form (^). In a Gaussian 
neighbourhood of F the metric h of S then has the form h = dx''^ + y^\h\g(x'). Now 

we introduce a scale factor a = 27r ^ Jo\h{z)\~^/'^ dz^ , where L denotes the length of a 
(projected) geodesic along all of orthogonal to the orbits F. Define a new coordinate 
X by x{x') = a Jq \h{z)\^^^^ dz, then we get a convenient representation of h as 

(9) h{x) = A{xf {dx^ + a^g{x)) , 

with A{x) = a~^|/i(a;)|^/^ defined on 5^, and g transforms under some element of GL(2, Z) 
like g in (^) . By the way, from this explicit formula and the definition of g in (||) it follows 
easily that the relation r = Aa holds. 

The Laplacian of S acting on functions ^ then has the form 

= -A-2(^" + A-i^V) + , 
with AV' = -h^^dAdsi^ since f% = 0. 

Finally we cast the second fundamental form k of S into a convenient form. Set K := 

k{m,m) = A~'^kii and observe kiB = — "Vbtt-i = —Ar]B we get 

(10) k{x) = A^K dx -2Ar]Bdxdy^ + KABdy"^ dy^ , 

where all quantities on the right-hand side depend on x. Taking the trace yields the 
relation 

H-K = tiK . 

2.1.3 The 4 geometry 

We will describe the 4-geometry (locally) in terms of a PMC foliation by symmetric 
Cauchy surfaces (S, h). It turns out, that there is not much left to do. Of course, we need 
the strong energy condition in M, and we can choose S without loss of generality to have 
non- vanishing second fundamental form. In view of theorem 2.2 of part I this is enough 
to guarantee the existence of a local in time PMC foliation of a neighbourhood of S in 
M. The remaining question is, if the leaves of the foliation turn out to be symmetric. But 



8 



the arguments given in the corresponding place in |H3| do apply to the present situation, 
since the Laplacian of S splits into the Laplacian of F and a part depending only on x, 
which coincides with the one in the plane symmetric case and we end up with 

2.1. Proposition. 

Let {M,g) be a globally hyperbolic, spatially compact spacetime with local U{1) x U{1) 
symmetry, obeying the strong energy condition. 

Then there exists a local in time PMC foliation {St} in M , covering a neighbourhood 
]ti,t2[x5^ of^ = So, in M. Moreover, ifT, is symmetric, then all the leaves of the 
foliation are symmetric, too, and there are coordinates (x^) = (t,x,y^,y^) adapted 
to the foliation, which cast the metric into the form 

1 2 A^u^ U2 vz\ 

^2 

g22 923 

Z^3 532 533/ 

,A _ jyA^ _ g^^D^ ^ g = r'^g are functions on ]ti,t2[xS^ , with 



(11) 
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v 



where N, A, v 



9ABV^ , 9 



r = Aa and a = a{t) only. The quantities with an overbar transform under some 
representation of GL(2,Z) after each transition x 1 — > x + 2tt in S^ . The shift 
functions are fixed by the conditions u^{t,0) = 27r) = and 0) = 0. 



The orthonormal frames and {V'^} associated with the 3 + 1-split read explicitly 

n = N~\dt-u'd,) 
m = A^"^ di 

62 = (Aa)-^ (+e-^/2pQgj^y/2 82 - e^l'^ ^xu^iV 83 

gij^i^ y/2 d2 + pQgj^ y/2 03 



eo 



62 



63 



V2 



63 = {Aay 
Ndt 

A{i^^ dt + dx) 

{Aa) pogh y/2 (p2 ^ ^^2^ ^ ^^^^ y^2 {u'^ dt + dy^] 

{Aa) gij^j^ y/2 (p2 ^ ^^2) ^ ^^gj^ y/2 (p3 dt + dy^] 



Finally, the Ricci rotation coefficients will be calculated. Some expressions are not ex- 
plicitly given, since they turn out to be too complicated. Instead, the dependence on the 
geometric quantities involved will be specified in square brackets: 

= A-^N-'^N' 



-K 
{Aa) 



*7Si 


= Ho 








= Ho - 


4^,0 




7ii 


= Hi - 




4^,B 


7lB 


= 7io = 




4^,B 


Vab 


= 7ao 




= 






Hb 


= Hb[ 






iQB 


= 'iU 


hj 





-1 — c 

eB Vc 



-CA^eB^ KCD 



'1 ,,1 



7oi 



where the dependence on the derivatives of V and W is linear. 
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2.2 The field equations 

Again we represent the matter quantities hy p = T(n,n), j = —T(n,m) = A~^ji and 
Sij — T'( Bj, Cj ). 

2.2.1 2+1 decomposition of the constraints 

First consider the Hamiltonian constraint 

iJ + F^ - = levrp . 

Decomposing R by the Gauss-Codazzi formula and k by ([l^) we get the following expres- 
sion: 

(12) m(trA) = 8^p--H'trK + |7/|2 + i(|(trA)2 + |Ap) + l(|(trK)2 + |K|2) . 
The first component of the momentum constraint 

V'kii - ViH = Snji 
can be decomposed into 

(13) m(tr k) = -Snji - F tr A + | tr k tr A + X^^kab 

(where the symmetry forced div rj = 0). The other components of the momentum con- 
straints are calculated analogously to (again some terms cancel out due to symmetry) 

(14) m(r?B) = -SirjB + (tr A)r?B . 

As in |H3| a convenient form of the equation is achieved by introducing the null expansions 



(15) := -(trAibtrK) = -A:±(r) 

r 



(compare the formulas for A and k in ^.l.lD , which yields an easy formula for the Hawking 



mass 



(16) M := miy(F) = -ir(ri?+)(ri?_) = -ir(^V"r)(*V„r) , 
by = 4/r2(''V"r)(^V«r). 

Adding and subtracting the constraint equations now yields 

(17) m{^±) = -87r{pTji)~\r]\'--'^\X±ii\^TH^±-^^l 

and the transition to the variables u;± := r'd± casts this equation into 

m{ijj±) = m(r)'d± + rm{'d±) 

= -87rr(p + ji) - r|?7|2 - lr|A ± ^ Huj± + ^ {uj+uj_ - 2^^ , 

where the identity m{r)'d± = ^((a;+ +ll!^)u>±) has been used. 
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2.2.2 2+2 decomposition of the field equations 

With the 2+2-geometry in mind, where we denote tensor components with respect to the 
frame of projection operators {qa,PA}, we write the equation ( [7a| ) for '^Rab in terms of g, 

'Rob = Rab - 2r-^YaYbr + \{YagAB){Yb~g'^'^) + 2g_ab\r]\^ , 

as a field equation in the quotient spacetime {B,g) for the radius function r. We stiU 
have to eliminate the unknown Rab = —Kgab, where K := -Roioi denotes the Gaussian 
curvature of {B,g) in {M,g). Contracting the equation for Rab and replacing Ar by (^) 
yields an equation for K: 

K = + 'Ri) + r-2(V'=r)(Ver) + \{Y'Ub){Yc~9'''') + 3|r?P 

Inserting this into the equation for *^Rab results in 

YaYbT = - \r-\Y''r){Ycr)9_ah - \r{^'Rab - \{'R-c - 'Ri)9_ab) 

+ \r{{Ya~gAB){Yb~9^'') - |(Y^9AB)(V,g^^)5a6) - hr\^?g_ab 

Now we are going to interpret the first three terms appearing on the right-hand side of 
the field equation: Recalling the definition of the mass function M in ( p!^ ) we can write the 
first term as y^gab- The curvature expression in the second term can be rewritten using 
Einstein's equations in terms of the projected energy momentum tensor T in the form 
^T^iXab - trTgab). 

To simplify the third term, there is more structure involved than already presented. The 
space of two-dimensional metrics Bil(2,R) is a three dimensional real vector space, whose 
elements we denote hy (p = (jf Ei, where the basis {Ei\ has been chosen as Ei = (g i )> 
E2 = ( -1 ) ) = ( 5 Q ). In this picture the determinant is a quadratic form on Bil(2, M), 
and the matrix representation with respect to the given basis turns out to be — "^/ij, where 
r]ij = diag( — denotes the standard Minkowski metric in three dimensional Minkowski 
space. The restriction of this space to the hypersurface r/((/), 0) = — 1 <^=^ detc/) = 1 is 
the well known hyperbolic plane H^^ and we have established a correspondence between 
the two-dimensional unit-determinant metrics g and the elements of the hyperbolic plane. 
Considering as a function on the quotient manifold -B, introduced in |2.1.1) we have a 
map 4> '■ iB,g) — > {H'^ C M?'^,7]ij), (t, x) i — > 4)^{g{t,x))Ei. Now let us interpret for a 
moment as a wave map from B to H^, then the energy momentum tensor '^T associated 
with <j) reads 

(19) ^Tab = Ya4>'Yb4>i - h{Y"¥Yc4>i)9_ab 

and with the correspondences (pi = r]ij(j)^ = —g^^ and (picf)^ = — ^ gAB g^^ , we see that 
the third term in the field equation for r can be rewritten as —\r'^Tab- 

Collecting all these things together we get the field equation for r in a convenient form: 

(20) YaYbT = ^9ab - 47rr (r,fe - tr Tgab) - \r ''Tab - \r\r]\^g_ab ■ 
Differentiating the mass function we find for the mass flux equation the formula 

(21) YaM = [Airr^iTab - tr Tgab) + jr ^Tab) (VV) + yH'^YaV . 
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The contribution of '*'T to the right-hand sides of these equations is known in so far, that 
"^T obeys the dominant and strong energy conditions. Furthermore, we can again rewrite 
its components using the formulas for the trace free parts of A and k in (remember, 
that Yo, Vi act as n, respectively m, on g, in the formula for "^Tab)- The result is 

(22a) V = V= |A|^ + |k|2 
(22b) ^ii = -2~XABfi^^ 



Now we want to analyse the underlying wave map. To make things more explicit, con- 
sider the parameters V,W of g as coordinates of in view of the correspondence 
^'£^j. The pullback under this parametrization yields a representation h of the metric 
iriij)\H^, with the explicit form hjj = dV"^ + (cosh F)^ dVF^, leading to 1^23 = tanhF and 
rig = — sinh V cosh V as non- vanishing Christoffel symbols. The pullback of the map cf) un- 
der this parametrization induces a map (j) '■ {B, g) — > {H"^, h), (j){t, x) = {V{t, x), W{t, x)). 
The wave map (j) is then defined to obey 

(23) V^Vc'/'^ + ^uY^^^Yc^'^ = . 



Let us turn back to the field equations. Writing equation ( ff^ ) in the form 

m-AB + 2rjAf}B = -\r-^gADY"{r^9^''YcgBc) = -l9ADY"ifYcr^9Bc) 
we arrive after some calculation at 

Y^ir^YcMB) =r^f{Y'~gAD){Yc~9Bc) 

- 167r(rAB - \ iTfgAB) - ^VAVB - ^Ivl'^gAB) , 

where T denotes the orthogonal projection of the energy momentum tensor T into {F,g). 
More explicitly we get 



(24a) 
(24b) 



or 



Y^ir^YcV) =r2 sinh V cosh V{Y^W) (YcW) 

- 2r2(coshy)-i((r23 - ^g^'^TABhs) - Umv3 - ^g^'^VAiiBh^)) 

Y^ir^YcW) = - tanhV{Y^W){YcV) 

- r\coshVr'[ie-'^T,, - e^Tsa) - U^-"^ iV2? - e^(%)^)) 



Y^YcV - sinh y cosh T/(V'=T^)(VcTy) = -2/r{Y^r){YcV) 

- 2(coshF)-i((T23 - y^'^TABhs) - h(.mV3 - y^^'vAVBhs 
Y^YcW + tanhV{Y^W){YcV) = -2/r{Y^r){YcW) 

- (coshy)-i((e-^T22 - e^Tas) - U^""^ {V2? - e^(%)^) 



(25a) 
(25b) 



The left-hand sides of equations (^) coincide with the the left-hand side of ( [2^ ) and we 
conclude, that V and W solve an inhomogeneous wave map. 
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2.2.3 3+1 decomposition of the field equations 



Now we want to calculate the ADM equations and bring them into a form most similar 
to the equations in the plane symmetric case of part I. The formulas concerning the 3- 



geometry in |2.1.2| are of particular importance for the calculations here, as well as the 



definitions of A, k and rj in |2.1.1| . A straightforward calculation then yields the following 
set of equations. 

The constraint equations 

(26) = 1^5/2 (^H^ - l(H - Kf - _ I ~ |2 _ |)^|2 _ 2|^|2 _ ig^^^ 

(27) K' = -3A~^A'K + A'^A'H + H' - Akab'^^^ + Svr^i 

(28) 7]^ = -2A~^A't^b - Svr^iB • 

The equations fixing the foliation 

2 



(29) N" = -A'^A'N' + A^N[\{H - Kf + K'^ + \k\'^ + 2|r/p + 47r(p + tr 5) ) - A 

(30) H = l + v^H' 

(31) {dt-v^d^)H' = u^'h' . 

The evolution equations 

(32) A = -NAK + Av^' + A'v^ 

(33) h = -'^Na{H -m)-av^' 
k =u^K' - A''^{N" - A'^A'N') 

(34) + 7V( - 4^-^/2(^1/2)// ^ (^-2^/)2 ^ 

- |A|2 + 2|?7|2 - 87r(A-25n + ^{p - tr S)) 



fjB ={\N{H + K)- 2A'^A'u^)i-iB - 2NkBcrf 



(35) + \A'^A'{H - K)uB - [k'ilDC + \{H - K)\bc - A~^A'hBc) 

- SttAi^^Jb 



As in |H3] integration of ( |33D over the circle yields the analogous equation for the first 
component of the shift vector, 



(36) v^' = -\N{H -2,K) + \ f N{H-m) 



and the definition of the second fundamental form provides some additional equations, 
namely 



(37) = -2NAr]^ 

and the coordinate representations of A and k in |T|. The remaining equations for V 
and W are still missing. To this end we supplement the system (p^)-([37|) by the field 



equations (24) or (|25|), and end up with the full system of equations. 

As in the plane symmetric case, we can formulate an existence and uniqueness theorem 
for solutions of the equations (^)-(|37|), (p5|). First we define a symmetric initial data set 
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for a spacetime with local U{l)xU (1) symmetry by the smooth collection (S, h, k), where 
S denotes a (possibly non-trivial) torus bundle over the circle and the fundamental forms 
h and k are represented in a suitable coordinate system as @ and ([lO|), respectively. If 
there are matter fields, then we assume the matter data and equations to be smooth and 
symmetric, leading to a well posed Cauchy problem coupled to the reduced field equations 
in harmonic coordinates. 

On the universal cover S of S the induced data then is invariant under the action of 
G, hence G acts isometrically on the whole Cauchy development, which induces a local 
U{1) X U{1) symmetry on the Cauchy development of (T,,h,k). Assuming A := + 
4:11 {p + tr S) > somewhere on E, we get for £ a unique symmetric local in time 
PMC foliation, defined on some time interval containing to, with E = StQ, and coordinates 



described in proposition 2.1 and we just have proved: 



2.2. Proposition. 

Let (S, h, k) be a symmetric initial data set for a spacetime with local U{1) x U{1) 
symmetry, with matter obeying the strong energy condition and A > somewhere on 
S. Further, let tQ denote an arbitrary real number. 

Then there exists a 6 > and a PMC foliated surface symmetric spacetime (M, g) 
dijjeomorphic to ]tQ — 6,tQ + dlxY^ with an embedding l : T, — > M, satisfying 
= StQ and i^^h, i^k are the first and second fundamental form of St^ in {M,g). 
(M, g) obeys the strong energy condition and g can be written in the form described 
in proposition WI. This construction is unique up to the choice of to and 6. 



2.2.4 The expanding model 

Now we can proceed in close analogy to the corresponding analysis in part I. The formula 
(1T6| ) for the Hawking mass shows, that gradr is timelike as long as M > 0. Indeed, 
proposition 3.1 in ||^] proves this, provided the spacetime is not flat and the dominant 
energy condition holds. 

Therefore, under these conditions gradr is timelike and ?9_ have fixed and opposite 
signs. Without loss of generality we choose the time orientation, such that gradr is 
past pointing and dr is future pointing (by the induced time orientation of the cotangent 
bundle). Then ■d^ > and < which classifies the spacetime as expanding in the 
sense described in part I. 

Again we expect the singularity in the distant past from our symmetric initial data Cauchy 
surface S and any symmetric Cauchy surface S" in M is not maximal with mean curvature 
not everywhere positive on S. 



2.3 A priori estimates for the field equations 

Assume the dominant and strong energy condition to be fulfilled in {M,g). Let S be a 
symmetric Cauchy surface in M with strictly negative mean curvature H and denote by 
{St}, t £]ti,t2[ the local in time PMC foliation with T, = Sq and the time orientation 
chosen in correspondence with |2.2.4 , such that H decreases with decreasing PMC time. 
We consider the past D~{T,) of S. In D^{T,) the mean curvature is bounded from above 
by some H < 0, and H = H only on S. Thus \H\ is bounded from below and as long as 
H remains finite, we find the following estimates. 
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Consider first the constraint equation in the form (^). The dominant energy condition 
gives the inequahty 



mi 



from which we infer the basic estimate 

(38) \^±\<4\H\ =^ \A-^A'\<C, \K\<C . 

Now we can perform nearly the same estimates, as we have done in the corresponding 



place in |H3|| . The additional terms appearing here in the more general equations (26)-(27) 



(p9|)~(p3D and (36) do not cause any problems. There is only one additional estimate, which 
will become important in our further analysis, resulting from the integration of equation 
(|2|) along 5"^ bounding not only f^^p, but also /jidAp + |kP + jr/p) = f^^[-t'p-\- |r/|2), 



and we can state the analogous 
2.3. Proposition. 

Let (M, g) be a globally hyperbolic, spatially compact spacetime with local U{1) x U{1) 
symmetry, obeying the dominant and strong energy condition. Assume the existence 
of a symmetric Cauchy surface E with strictly negative mean curvature. In particular 
we get from proposition \2. 1\ a PMC time coordinate t, ranging in ]ti,t2[ with S = 
{t = 0} and H decreases with decreasing t. 
Then we have uniformly on ]ti,0] 

|A|,|A-^|,|A'|,|a|,|a~^|,|i/|,|if'|,|E:|,|iV|,|A^-^|,|iV'|,|i/^|,|i^^'| < C . 



Until now the control over the coefficients of the fundamental forms h and k is not complete. 
We still need bounds for the metric coefficients V, W, the components rj and k oi k and 
the remaining components of the shift vector. We will get bounds for most of these 
quantities here in this section, using the bound for fgi{p + "^p + |??P)- 

2.4. Proposition. 



Under the hypotheses of proposition 2.S we get uniform bounds for 
|F|,|l^|,|ryB|,k^|>^'| 

on 0]. 



Proof. 

The bounds for 77^, and are simple consequences of the bounds for V and W , 
since having bounded them, we can conclude as follows: The bounds for V and W allow 
coordinate components to be bounded in terms of components according to the orthonor- 
mal frame vectors 62,63 in {F,g) defined in We adopt the convention, that a hat 

above indices denotes components with respect to an orthonormal frame. Having this 
in mind we see, that J^i p bounds J^i js < C J^i by the dominant energy condition, 
so integration of the constraint ( p8| ) yields a bound for the difference A^tjbI^I, '^^lich is 
independent of t. Using tjb < Crjj^ < C(l + |?7p) we can bound the integral fgi\A'^r]B\ by 
C J^i(l + |r?P), which together with the estimated difference yields the desired estimate 
for r]B. Now we get immediately a bound for u^' by inspection of equation (|37|), and the 
condition i'^{t,0) = bounds u^. 
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1^-2 (y'2 + (cOshF)2^'' 



Let us now investigate the field equations (24) for V and W, considered as equations on 
{B,g) endowed with the coordinates {t,x), as described in the end of |2.1.1| . Then we 
calculate explicitly in these coordinates 

|A|2 + \k\' = i(m(5AB)m(5^^) + nCgAsMg^'')) 
_ 1 

2 

+ ^N-^ i^{V - vWf + {coshV f{W - v^W'f 

and we see, that the first term on the right-hand side of equations ( p4| ) can be bounded 
in terms of this expression. 

For the second term on the right-hand side of ( p^ ) we exploit the special structure ap- 
pearing there. The first aim is, to express the tensor components with respect to an 
orthonormal frame. Looking at the definitions of the coframe {ct^} defined in |2.1.1| we 
recognize the relations 

e-^/'^dy^ = r-^ (- sinh V/2 + cosh V/2 a'^) . 

If we mark tensor components expressed in this basis by a tilde, we see that the terms in 
brackets in both equations ( |2^ ) are easily rewritten with respect to this basis, by simply 
putting a tilde above each tensor component and deleting the factors e^^^'^. But we have 
schematically r/ = rjBdy^ = fjE e^^^'^ dy^ = fB{cosh.V/2)fjB o'^ , with /^(cosh l//2)f/B = 
77^ is the component according to the orthonormal frame. A similar relation holds for 
Tab- Since we have rj^ < 1 + |r/|2 and < p we conclude, that the terms in brackets 
in ( p4|) are bounded by C cosh V{1 + p + which leads together with our result about 
the first term on the right-hand side of (|2j) to 

I V^(r2v,0)| < C(l + p + I A|2 + \k\^ + |r/p) , 

and since the integration of the right-hand side along is already bounded we get 

/ |V^(r2Vec/>)| < C 



uniformly for each t G 0]. Note, that we use the wave map (p = (V, W) defined in 2.2.2 
in order to abbreviate some formulas only. 

On the other hand, we can interpret the integrated quantity as a divergence, a fact we will 
take advantage of to get rid of the integral sign. To this end consider a point (t, x) in the 
quotient manifold B, with t G]ti,0[. We need estimates uniformly in t for t approaching 
ti. We define the (upside down) characteristic triangle T in B, by its counterclockwise 
oriented boundary dT = 7+ -|- 70 + 7_, where 7-1- denote the characteristic curves of (24), 
with 'y± = k± = m±n connecting (i, x) with (0, x±) and 70 is the curve in the hypersurface 
{t = 0} from (0, to (0,a;_) with 70 = — m. Define v £ T^F by v"" = r^V'^c/), then the 
corresponding 1-form oo reads w = i{v)^, where VL denotes the volume form with respect to 
g. Expressing tensor components with respect to the projected orthonormal frame {aa = 

qa} (see 2.1.1 ) we have v"" = r"^ i )' ^"^^ ^ ~ {rn{4>)q^ + n{(f))q}^ . The pullback 



of uj along dT has the three parts 7^.0; = k±{4> o j±{u))du, ^qUJ = r n{<j){x))Adx. 
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Thus, setting r := -J2t 



UJ 



IdT 



r'^k+{(j){j+{u)))du+ / n{(f){x))Adx + / +r^A;_(0(7_(n))) 



r'^k+{(l){j+{u)))du 



-r^(f) 



7+(0) 
7+(t) 



r 



7-{0) 
7_(t) 



'0 



+ 



+ 



n{4>{x))A dx 
n((/)(x))yl dx 



where integration by parts has been carried out in the last step. Note, that the last term 
is bounded, since that integration takes place on the line t = in i?, where (j) is smooth. 
In addition fe±(r^) is bounded everywhere by proposition p.3| . Further, the first terms can 
be evaluated to — — ?'^</'|!^~(°j = 2(r^</))(t, x) — (r^(/))(0, x+) — (r^0)(O, where 
the last two terms live on the line t = 0, and therefore are bounded. 
Stokes' theorem now applies to the present situation, jgrpOJ = jrpdw = Jrp{d\vv)Q-. 
Putting all this together we get an estimate 



\rmt)\\i 



< C 1 + 



is)\\L^ ds + 



where the last term is already known to be bounded uniformly in t by our previous 
analysis. Hence we can apply Gronwall's inequality and the bounds for r, obtained 
in proposition 2.3 complete our argument, arriving finally at |1/|,|VF| < C uniformly on 
]ti,0[. ■ 



Unfortunately, there is still a lack of control over the components k k (respectively its 
traceless part, since (tr k) g = (H — K)g < C). It turns out, that the necessary bounds 
for the first derivatives of </> depend on a bound for the matter quantities and vice versa. 
The analysis of this dependence will be performed in the next section, where we have to 
take a specific matter model into account. 



2.4 Higher order estimates 
2.4.1 First order estimates 

We start with the derivation of first order estimates for the metric coefficients = {V, W). 
To keep things simple, we consider collisionless matter, as described by the Einstein Vlasov 
system. 

The particle density / on the mass shell is governed by the geodesic spray X with 
X = v^'^Ca — '^I'^u'^^'^'^ the killing vector fields Sa define two conserved quan- 

tities 5f(X, 5^), which bound v"^. To bound we consider the characteristics of the 
Vlasov equation (compare the explicit formulas in ]H3|] ): 

<C{l + Qi{s)+v\s)) , 

where Qi := D := {dt,dx) and we used the fact, that in view of the boundedness 

of v^, the vanishing of *7q]^ and the special form of the non- vanishing rotation coefficients 
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no term involving the product Qiv"^ occurs. Thus we find for the quantity Pf{s) = 
{sup|u| I V G supp/} (which measures the matter quantities, compare |H3| ]) the Gronwall- 
hke estimate 



1 + Pf{t) <c(^l + Pf{0) + j\pM + Qi{s)) ds^ 



This inequality shows, that we indeed have to estimate the matter quantities together 
with the second fundamental form k. The next step consists in finding a complementary 
inequality for Qi, which yields in combination with the inequality just obtained a true 
Gronwall estimate for 1 + Pf + Qi. 

We need again the field equations for V and W, but now it is more convenient to analyse 
the equations in the 'wave map form' (^). First consider the left-hand side of the field 
equation. We want to express this derivative operator in terms of the characteristic vector 
fields k± = m ziz n tangent to the characteristic curves 7± introduced in the previous 
section (in the proof of proposition In the quotient manifold {B,g) one calculates 

the non-vanishing Ricci rotation coefficients to 

foi=fo^ = {AN)-^N' 

f,,=f,, = {AN)-\A-{Av')') , 

thus Yea§.b is bounded. The characteristics k± expressed in the coordinates of B are simple 
linear combinations of the Ca, hence Yk±k± is also bounded in B. We will use this fact 
to estimate commutators such as [m, n] or [k^,k^]. Transforming from the orthonormal 
frame {sa} in {B,g) to the frame {k±}, we can write the left-hand side of (p5D as 



V^V,0^ + f fjVy Ve</>^ = Vfc^ (A:-(0^)) + rfj k+{(l)')k_{<p') + [k+, k^icj)'') 
and there is an analogous equation for kj^ and A;_ interchanged. 

Let us turn now to the right-hand side of (p5|). The terms not known to be bounded are 
of the form k±{(f)) and Tab- Fortunately we have control over the matter term, since the 
bounds on V, W and v"^ justify the inequality Tab ^ C Pf and we obtain 

Yk^ {k.icp'')) + tfjk+icj)')k.{cp') < C(l + A:±(<A) + Pf 

The left-hand side is still non-linear in the first derivatives of <j), so this inequality is not 
in the appropriate form to apply some kind of Gronwall estimate. But we can overcome 
this difficulty by adapting an observation of Gu ( |Gu| , see also |^ for explanation), who 
considered wave maps defined on two dimensional Minkowski space. We already performed 
the first step, consisting of a transformation to characteristic coordinates. Now we define 
the vector field k± over the map (p, 

C\ C) C) 

k± := Mk±) = k±icP) = A:±((^^)— ^ = k±{V)— + k±{W)- 



thus bounding the length of k± accomplishes the estimate of Qi. The vector field k± over 
4> \s a. section in the bundle (ffiTH"^) over B and the covariant derivative operator (j)*V 
over (j) acts like {cj)*V).uX = v{w^)-^^ -\- 'W^'^,f),v-Q^ for every vector v £ TB and every 

X = w^-^ with : B ^ H^. 

Evaluating this expression with v = k^ and X = k- gives 

(rV),^fc- = k+ik^{<P^))^ + tfjk+i^^)k.icpJ)^ 
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(and analogously for k± interchanged). The components in this equation are similar to the 
left-hand side of the inequality already obtained for the field equations. The remaining 
term involves k±{4>) times the connection coefficients in B with respect to the frame {k±}, 
which we know to be bounded. Thus we get for the field equations the inequality 

\{cl)*V)k^k4<c(l + k±{<p)+Pf) , 

expressing an estimate about the growth of k± during the transport along the characteristic 
curve 7^. This inequality is in the appropriate form for a Gronwall-like estimate: Taking 
the maximum for each fixed t allows us, to combine the inequalities for and |A;-|, 
replaced collectively in terms of Qi. Hence 

Qi{t) <c[l + Qi(0) + j\pM + Qi{s)) ds 

where 7-1- has been reparametrized in terms of t. Combining this with the inequality for 
Pf we arrive at 

1 + Pf{t) + Qi{t) <c(^ + Pf{0) + Qi(0) + _^ (1 + Pfi^) + Qii^)) ds 

Performing a Gronwall argument we have proven 
2.5. Proposition. 

Let {M, g, f) be a globally hyperbolic, spatially compact solution of the Einstein- 
Vlasov system with local U{1) x U{1) symmetry, which possesses a symmetric Cauchy 
surface S with strictly negative mean curvature. The PMC time coordinate t ranges 
in ]ti,t2[ with T, = {t = 0} and H decreases with decreasing t. 
Then we get uniform bounds for 

\V\,\Wl\V'\,\W'\,p 

on 0]. 
2.4.2 Second order estimates 

We still have to do one further step before it will be possible to apply some iteration 
scheme. We need to establish bounds for the second derivatives oi (j) = (y, W) together 
with bounds for the first derivatives of the matter variables. Again it turns out, that it is 
not possible to get separate estimates for these quantities. 

First we differentiate the field equation (|2^) and obtain an equation of the form 

k± {k^id^cj)^)) + 2tfj k±{(j)^)k^{d^(j)'^) = [...]d^TAB + lower order terms , 

where [...] is an abbreviation for some term involving already bounded quantities, in par- 
ticular we see by a quick look at the 3+1-field equations, whose right-hand sides are now 
bounded, that for the Ricci rotation coefficients dxj < C holds. Note, that the differenti- 
ation kills the nonlinearity in the equation, but we have to deal with first order derivatives 
of the matter variables instead. These quantities cause serious trouble, because if we 
differentiate the Vlasov equation, terms involving the second derivatives of (j) times first 
derivatives of / come up, thus there is no direct Gronwall argument possible. To attack 
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this difficulty we have to combine the equations, using an idea of [GS|. 



First we integrate the equation along the characteristics 7-1- to get an integral inequal- 
ity. In order to apply a Gronwall argument all what remains to do is to bound the term 
Xy^ QxTab- Then we express Tab in terms of its components with respect to the orthonor- 
mal frame, producing merely some already bounded quantities and the frame components 
of the energy-momentum tensor for Vlasov-type matter look like J^^ (J va^b/v^ dxf dv) . 

Finally we represent dx by a linear combination of k± and X := dt + {N A^^v"^ /v^ — i^^)dx, 
which is the part of the characteristic vector field for the Vlasov equation, lying in B. 
The transformation to this basis is obviously bounded and we can proceed as follows. 
The part involving k± permits a direct application of the integration by parts rule, which 



contributes something bounded by proposition 2.5. The remaining part can be treated in 



a similar fashion, after inserting the Vlasov equation into X{f). This yields terms involv- 
ing only first order derivations of (j) (arising from the ^7's) times and again we can 
perfom an integration by parts with respect to the velocity integral. The result consists 
in terms already bounded and we are through. Applying now Gronwall's inequality we 
arrive at a bound for ||93;-D</>||2,oo and inserting this result into the field equation (25) we 
get Q2 ■= \\D'^4'\\l°° bounded. 

An immediate consequence is, that we also have established a bound for the first deriva- 
tives of /. Differentiating the Vlasov equation with respect to x or yields an equation 
for dxf or d^f respectively, with bounded characteristics and an inhomogeneous term, 
consisting of (^e^i)", ''7^1^ and their derivatives with respect to x. These terms are either 
bounded by what has been said above or by inspection of the field equations (p6D-(|37[). 
Having bounded the spatial and velocity derivatives of /, the structure of the Vlasov 
equation bounds immediately dtf and we have proven 

2.6. Proposition. 



Under the hypotheses of proposition 2.5 we get uniform hounds for 
\V\,\V'l\V''\,\W\,\W'\,\W''\ and \p\,\p\ 

on 0]. 



2.4.3 The iteration scheme 

Following the analysis of part I we need some additional framework. 
2.7. Definition. 

T:= {A,a,V,W,N,i^\i^^,H,K,7]B) 

It turns out, that what we have done in the previous subsections are the first steps towards 
the matter regularity property of part I. Here we will define this property adapted to the 
present situation with the notational conventions introduced as follows: a denotes a multi 
index for derivatives in {B,g) and D := {dt,dx)- Then the matter regularity property 
reads 

(39) < C =^ \D"^\ < C 

and we state 
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2.8. Proposition. 

Under the hypotheses of proposition \2.3[ ( |39D holds for all a uniformly on ]ti,0]. 



Proof. 

We prove the statement by induction with respect to \a\. The induction hypotheses 



for a = is contained in the statements of the propositions 2.5 and 2.t. To proceed 
further let us assume \D°'J^\ < C for some a with \a\ =: p > 1. Then we can assume 
< C by induction. We have to show l^^'I'l < C, or extending our previous 
notation, Pp := ||I?p/||loo < C and Qp+i := < C. 

Fohowing the analysis of the previous subsection we set q := p — 1 and differentiate the 
field equations ( p5|) with dxD'^, which yields schematically 

k± {k^{d.,D^)) =[a,^«7, d^DPr, d^Dm {k^O^D^) 

+ [D^] d^D'^TAB + [d^D^, d^D'^r], W^Tab] , 

where quantities in square brackets abbreviates some expressions formed by them, whose 
detailed structure is not important for our analysis here. We want to perform the same 
kind of argument as we have done in the proof of the second order estimates. In order to 
do this we must bound the quantities in the square brackets: 

• The terms in the second and third square bracket are bounded by the induction 
hypotheses. 

• For the first square bracket we can proceed as follows. is already bounded. 
For dxD^r we have to estimate dxD^A = d^D'^A. This can be done by applying D*^ 
to ( p6[ ) in a strictly analogous manner as in the first item in the proof of lemma 3.6 
in pi3| , applied to equation (p^). The definition of 7 in p. 4.1 shows, that the terms 



in dxD'^^ not already bounded are dlD'^f^ , d'^D'^N and dxdtD'^A. Applying to 
(^) gives a bound for the first quantity. Using D'^ on (|29|), then the argument in 
the fourth step in the proof of lemma 3.6 in p^3| |, applied to ( p9[ ) bounds the second 
quantity. Now all quantities appearing on the right-hand side of {dxD'^ applied to) 
(^) are bounded, and this bounds dtdxD'^A. 

Turning now to the Vlasov equation and applying dxD'^ yields an inhomogeneous equation 
for dxD^ f with the same characteristics. Thus we can apply the same trick as in the 
previous subsection, substituting dx by k± and X. Then again we are concerned with 
integration by parts, which yields only bounded terms, and the inhomogeneous term of 
the differentiated Vlasov equation, which is also bounded by induction hypotheses. All 
together we can apply Gronwall's inequality to find H^x-D^c/iIIloo < C, and immediately 
Qp+i < C* by inserting the spatial bounds into the field equation. The Vlasov equation 
for dxD'^f now bounds dxD'^f and analogously for d^D'^f, which automatically bounds 
dtD^ f by inserting the spatial and velocity bounds into the differentiated Vlasov equation. 
Thus we have also bounded Pp, which completes the proof of the proposition. ■ 

To extend the local in time PMC foliation, our aim is to find uniform C°° bounds of all 



geometric and matter quantities, when the PMC time t approaches ti. Property (39) 
shows, that it is enough to bound the geometric quantities jy^T for all a as long as t or 
respectively the mean curvature remains finite . Therefore all we need is the analogue of 



the lemmas 3.6 and 3.7 in |H3], starting with the propositions 2^, 2^, 2^ and p.6| 



The first part is straightforward, given an arbitrary multi index a and \D°'T\ < C we 
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find \dxD'^T\ < C by inspection of the relevant (differentiated) field equations (|26|)- 
(compare the arguments given in the proof of lemma 3.6 in p3] where necessary), with 



the terms on the right-hand sides bounded by property (pQl). Moreover, proposition 
provides not only bounds for the spatial derivatives, but also for the time derivatives of (j) 
by the definition of "^/o. This allows us to proceed as follows: First we bound, starting from 
D°'J^ all spatial derivatives d^D"J^, k = 1,2,... Then we bound, starting successively 
from d^D'^J-, k = 1,2,... the quantities D^J^ for all multi indices /3 with = |a| + k, 
redistributing the spatial derivatives into time derivatives. This procedure will successively 
bound all derivatives of taking advantage of the fact, that in each step all lower order 
derivatives with the same order of time derivatives have been bounded as well as at least 
one more spatial order derivative in the step with one order less in time derivatives. In a 
more compact formulation, we have to show the boundedness of dfD^J^, having already 
bounds for dxD'^J' and D^T . 

We accomplish this step for each member of T ^ by again considering some of the (dif- 



ferentiated by L>") 3-|-l-field equations. We see immediately from (|3^ ) and (|33[) , that 
we have bounds for dtD'^A^ dtD'^a. The bounds for dtD°' applied on H, K and r/ are 
straightforward, too. Of course, dtD"(p is already bounded by property (|39D. Moreover 
we can strengthen the regularity. Since we have already a bound for dxD'^J' proposition 
|2.8| provides us with bounds for d^D^ applied to p and Dcj). Applying d^D" to the field 
equation (^) we first get a bound for dxD'^Dcf) and then for d^D'^cj) by inserting the first 
result into the field equation. This in turn used together with the bounds for the differ- 
entiated matter variable in the differentiated Vlasov equation, bounds dtD^p. 
Now we turn to the analysis of the differentiated lapse equation. First we see, that the 
bound for dtD^u follows from the bound of dtD'^N. For the latter one, we follow the 



argument given in the corresponding place in the proof of lemma 3.7 in |H3|, where all 
that is needed, has already been bounded by the arguments just given here and we are 
done. 

Therefore we end up with 
2.9. Theorem. 

Let (M, g, f) be a globally hyperbolic, spatially compact solution of the Einstein- 
Vlasov system with local U{l)xU{l) symmetry, which possesses a symmetric Cauchy 
surface T, with strictly negative mean curvature H < H < and H = H somewhere 
on S. 

Then all of the past of S admits a PMC foliation {St}, where t takes all values in 
the interval ]— oo,0] and H takes all values in ] — oo,^]. 



2.4.4 Improving the result 



Here we will try to do the same construction as in the corresponding place in |H3|, to get 



rid of the restriction concerning the fixed sign of the mean curvature H on the Cauchy 
surfaces. So we assume M to be non-flat and denote by S an arbitrary symmetric Cauchy 
surface in M. In order to follow the steps performed in the plane symmetric case in part 



I, remember first from |2.2.4| , that the spacetime is expanding. Note further, that the 
Einstein- Vlasov system fulfills the dominant and strong energy condition as well as the 
non-negative pressures condition. 

Then we find for the past domain of dependence D~{Y^) 
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Since the spacetime is expanding, dr is future pointing, thus r is bounded in D~{T,). 
By assumption, is bounded on the compact surface E. Unfortunately, unlike 
in the plane symmetric case we cannot conclude from the non-negative pressures 
condition and the fact that dr is future pointing, that dM is past pointing, since in 
equation (^) for the mass flux the term involving r] contributes with the wrong sign 
(luckily, the contribution of the energy momentum tensor ''T ([l^) for the wave map 
(/> does not cause any trouble, due to the energy conditions automatically fulfilled by 

To overcome this difficulty we impose the Gowdy-type symmetry condition, thus r] 
vanishes identically in the spacetime (and by the gauge condition z^^(t,0) = this 
is also true for P in view of (|37|) ). Then we conclude, that is bounded in Z)~(S) 
as desired. 

Now we want to apply theorem 2.1 in Q to establish a bound for the length of all 
timelike curves in D~(E). For this we have to adapt the proof a little bit. Inspection 
of the proof shows, that the argument relies on the inequality r < — ^ satisfied by 
the area radius along timelike geodesies, where the dot denotes differentiation along 
the curve. But looking at the field equation (|20| ) for the area radius easily establishes 
this relation (since rj vanishes) and we are through. 



The basic estimate ( |38D together with ( |16D shows, that r and M are bounded 



Lemma 2 in [BR] applies and we get bounds for the volume and its inverse for any 



Cauchy surface in D (S). 

The volume V(t) for the leaf St of the PMC fohation is given by V{t) = (47r)^a^^ J^i 
and a closer look at the estimates given in the corresponding place in [|H3| ] shows, that the 
arguments apply literally. This relies on the facts, that on the one hand the additional 



equations (25) are unaffected by the reparametrization of the foliation, thus the estimates 
for (p and its derivatives done in the previous subsections hold. On the other hand the 
construction in part I is based mainly on the structures introduced by the matter regularity 
property and the lemmas 3.6 and 3.7, a program we adapted successfully to the more 
general situation here. 

So we finally arrive at 

2.10. Theorem. 

Let (M, g, f) be a globally hyperbolic, spatially compact solution of the Einstein- 
Vlasov system with Gowdy-type local U{1) x U{1) symmetry and T, be a symmetric 
Cauchy surface. 

If (M, g) is non-flat then we can foliate all of the past of S by PMC hyper surf aces, 
where the time function takes on all values in the interval ] — oo,0] and the mean 
curvature of the leaves tends uniformly to —oo for t — oo. 

Using the PMC leaves as barrier surfaces we get the 

2.11. Corollary. 

In the situation of theorem 2.1(\ D^{T,) possesses a CMC Cauchy surface for each 



value of the mean curvature in ] — oo, mhiY. II[. 
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3 Conclusion and outlook 



We have seen in this second part, that the program initiated in part I has been successfully 
extended to the more general case of local U{1) x U{1) symmetric spacetimes in close 
analogy to the plane symmetric spacetimes. Therefore all what has been mentioned in the 
corresponding place in p^ ] applies. 

A corresponding analysis for spacetimes with local U{1) x U{1) symmetry has been inde- 
pendently performed in the work |Xn| of Andreasson. His construction leads to stronger 
results (as the work of Rein in the surface symmetric case of part I), but the time 
functions used are defined in terms of the symmetry. Thus it is not clear, how to general- 
ize them. In view of Andreasson's work, the present analysis can be seen as a suggestion 
pointing in a slightly different direction. 

In comparison with the spacetimes considered in part one the control of the momenta 
of the Vlasov particles turned out to be more complicated, in particular the coupling to 
gravitational waves required second order estimates before the iteration procedure could 
be performed. In this process we took advantage from the wave-map structure of the 
dynamical part of the geometry driven by the simple form of the Vlasov equation. 
Nevertheless one can hope, that the approach may be generalized to other matter models 
and geometries, since the obstructions seem to be of technical nature only. 



Acknowledgements: I want to thank Alan Rendall, who guided me through these parts 
of my work, which depend heavily on his previous developments. 
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